UNCLASSIFIED 


ad  4  2  5  131 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


NOTICE:  Wien  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  io  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


RESEARCH 

^  REPORT  219 

lO  NOVEMBER  1963 


IYPJMKS  OF  TOWED  UNDERWATER  VEHICLES 


Subprojtet  No.  SF  011  02  12  •  Task  No.  2563 

A.  G. STRANDHAGEN 
C.  F.  THOMAS 


UNITED  STATES  NAVY  MINE  DEFENSE  LABORATORY 
PANAMA  CITY,  FLORIDA 


I  ur  1 

i, 

U{» 


I 


COPY  HO.  51 


ABSTRACT 


This  report  presents  the  analysis  of  a  system  which  Includes  a  - 
maneuvering  ship  towing  an  underwater  vehicle  at  the  end  of  a  long 
flexible  cable.  The  equations  of  motion  for  both  the  cable  and  the 
underwater  vehicle  are  also  presented. 

The  cable  is  imagined  to  consist  of  many  interconnected  short  rigid 
segments.  The  equations  of  motion  for  the  system  are  formulated  twice 
on  the  basis  of  two  hypotheses:  first  for  a  simple  hypothesis  regarding 
the  inertia  of  an  accelerated  body  in  a  fluid,  and  secondly,  for  a  more 
complete  and  a  more  accurate  hydrodynamic  hypothesis. 
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LIST  OF  SYMBOLS 


Xx  or  V 


the  dot  over  a  quantity  denotes  differentiation 
with  respect  to  t,  i.e.,  d/dt. 


two  notations  are  used  to  denote  vectors:  either  the 
vector  bar  over  the  letter,  or  a  subscript  to  denote 
the  component.  For  example  x  ■  (x^,y^,z^)  "  • 


Of.gj Y,p  these  Greek  letters  have  specific  meanings  assigned 

further  on  in  this  section.  They  else  are  dummy 
variables.  The  quantity  p  is  also  the  mass  density 
in  Appendix  E. 


sin(vx  ,  vs) 


ax 


A* 

v+l 


the  sine  of  the  angle  betwee  •  the  vectors  of  vx  and  va  . 

v\\ 

a  computational  abbreviation  for  + 


X  *  1*2,  • • « *  v  . 
a  computational  abbreviation  for 


e^) ,  where 


w 


g  '  v+l  f  dv+5)  d 


dvfl  ‘  dv+2 


v+l 


)  +  <  Vl/2) 


aJf2 


A  computational  abbreviation  for 
+  d  )■  +  (e  /2) 

g  '°v+l  r  av+5;dyfl  ^  v+r  ; 


^i 


component  of  moment  on  the  X  th  element  in  the  1  th 
direction  due  to  the  acceleration  of  the  fluid. 


b 


X 


the  buoyant  force  per  unit  length  of  the  X  th  segment 
due  to  the  displaced  water. 


B 


Xi 


component  of  buoyant  force  on  the  X  th  segment  of  cable. 


BM 


1 


component  of  moment  about  the  c.g,  of  the  vehicle  due 
to  the  buoyant  force  acting  on  it. 


b 


Xij 


a  computational  abbreviation  for  e^s^s 


Xj 


Cj.  the  coefficients  of  the  linear  orthogonal  coordinate 

J  transformation  connecting  fixed  inertial  coordinates 

with  coordinates  fixed  in  the  towed  vehicle  with  their 
origin  at  its  c.g. 

q  the  proportionality  constant  between  the  square  of  the 

velocity  of  a  cable  segment  and  the  tangential  force 
per  unit  distance  due  to  the  velocity  of  the  water. 

C 

1, v+1  for  the  vehicle  is  the  same  quantity  as  is  for  a 

cable  segment. 

q  the  proportionality  constant  between  the  square  of 

the  velocity  of  a  cable  segment  and  the  normal  force 
per  unit  distance  on  the  segment  due  to  the  velocity 
of  the  water. 

C0  for  the  vehicle  is  the  same  quantity  as  C0  is  for 

a  cable  segment. 

C,  the  proportionality  constant  between  the  lift  on 

the  vehicle  due  to  the  fins  in  the  vertical  plane, 
i.e.  the  horizontal  deflecting  force,  and  the  square 
of  the  velocity  through  the  water  times  the  function 

FLVP(“  ‘  • 

C4  the  corresponding  constant  for  the  drag  on  the  fins 

in  the  vertical  plane. 

Cg  the  proportionality  constant  between  the  lift  due 

to  the  horizontal  fins  and  the  square  of  the  velocity 
through  the  water  times  Fj^O  “  3i). 

Cq  the  corresponding  constant  for  the  drag  of  the  fins 

in  the  horizontal  plane. 

Cl^  added  masses  i  «  '  j  £  3  (Reference  3). 

C2ij  added  masses  j+3>  i  +  3>3,  j  +  3>3  . 

03^  added  masses  ^  i  <  3,  j  +  3  >  3;  or 

^1+3, j’  1  +  3  >  3»  i  *  3  • 

c.b. ,  c.g,,  are  centers  of  buoyancy,  gravity,  mass. 

c.m. 
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v+i 


'ijk 


Xi 

* 

FXi 


A\i 


VX1 


GU 


(FIN), 


(FINM) . 


*Fapp\i 


<Facc>xi 


fdhf^  ^  ' 


is  the  length  of  the  X  th  segment,  for  X  *  1,2,...,  v  . 

distances  measured  on  the  vehicle  as  indicated  in 
Figure  A3,  where  i  «  1,  2 . 6. 

the  skew  symmetric  three  index  symbol. 

the  proportionality  constant  (Reference  1)  between  the 
normal  component  of  acceleration  of  a  cable  segment  and 
the  force  on  it  due  to  the  accelerated  motion  of  the 
fluid. 

equals  e  divided  by  the  mass  of  the  cable  element. 

the  total  force  applied  to  the  X  th  element. 

the  total  force  applied  to  the  X  th  element  exclusive 
of  the  mechanical  forces  at  either  end. 

the  sum  of  all  noninertial  forces  on  the  X  th  element 
which  depend  on  the  acceleration  of  the  element. 

the  sum  of  all  forces  on  the  X  th  element  exclusive  of 
those  on  the  fins  which  depend  on  the  velocity  of  the 
element  but  not  on  its  acceleration. 

the  sum  of  all  forces  on  all  the  X  th  element  which 
depend  on  neither  the  element's  acceleration  nor  velocity. 

the  force  exerted  by  the  fluid  on  the  vehicle  through 
the  vehicle  fins. 

moment  exerted  by  the  fluid  on  the  vehicle  through 
the  vehicle  fins. 

that  approximation  to  the  force  on  the  Xth  element  due 
to  the  velocity  of  the  fluid  obtained  by  neglecting  the 
rotation  of  the  X  th  eier-ent. 

the  force  on  the  X  th  element  due  to  the  velocity  of  the 
fluid  when  the  rotation  of  the  element  is  taken  into 
account. 

an  empirical  function  giving  the  lift  of  the  horizontal 
fins  when  the  vehicle  is  parallel  to  the  stream  whose 
velocity  is  unity. 
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flvf<0,“0*) 


4,m,n 


an  empirical  function  giving  the  lift  on  the  vertical 
fins  when  the  vehicle  is  parallel  to  a  stream  whose 
velocity  is  unity. 

gravitational  acceleration. 

a  computational  abbreviation  for 

ex  V'V^2’  +  W1  +  \ 

+  C^d^ia^  x  u^)q^  +  FIN  ,  where  FIN  is  the 
zero  vector  unless  X  *  v  +  1. 

is  the  angular  momentum  of  th  X  th  element  about  its 
c.g.,  where  X  ■  1,2,  ...,  v  +  1. 

equals  EtAnOx^x,,  where  (Am)  has  an  i  th  coordinate 
distance  x^  from  the  c.g.  of  the  element,  and  summation 
is  carried  out  for  all  Am  into  which  the  X  th  element 
is  decomposed. 

are  the  products  of  inertia  of  the  vehicle  about  its 
c.g.  referred  to  a  frame  of  reference  fixed  in  the 
vehic le . 

in  Appendix  E,  a  unit  vector  along  the  fixed  xx  axis. 

in  Appendix  E,  a  unit  vector  along  the  fixed  Xg  axis. 

in  Appendix  E,  a  unit  vector  along  the  fixed  X3  axis. 

in  Appendix  D,  a  proportionality  constant  which 

cancels  out. 

unit  vectors  along  ,  £e  ,  §3  axes  respectively, 
the  mass  of  the  X  th  element. 

the  sum  of  all  the  moments  acting  on  the  \  th  element, 
about  the  i  th  axis. 

the  sura  of  all  the  moments  on  the  X  th  element  about 
the  i  th  axis  which  depend  on  neither  the  element's 
velocity  nor  its  acceleration. 
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Vi 


the  sum  of  all  the  moments  on  the  X  th  element  about 
the  i  th  axis  which  depend  on  the  element's  velocity 
but  not  on  its  acceleration. 

the  sum  of  all  of  the  moments  on  the  X  th  element  which 
depend  in  any  way  on  the  acceleration  of  the  element. 

a  unit  vector  along  the  £3  axis  fixed  in  the  vehicle. 

a  unit  vector  along  the  intersection  of  the  plane  normal 
to  the  X  th  element  and  the  plane  containing  two  lines, 
one  along  the  axis  of  the  element  and  the  other  in  the 
direction  of  the  velocity  of  the  fluid  past  the  element; 
where  X  ■  l,2,..,,v+l. 

unit  vectors  related  to  the  direction  of  flow  of  fluid 
past  the  vehicle  fins,  where  X  »  v  v  2,  v  +  3,  v  +  4, 
v  +  5. 

a  general  displacement  vector  used  in  Appendix  E  to 
relate  the  Kirchhoff  force  to  an  inertial  frame. 

a  unit  vector  along  the  X  th  element,  X  *  1,2,..., 
v  +  1. 

,£or  the  towing  vehicle,  a  unit  vector  corresponding  to 


the  mechanical  force  exerted  on  the  X  th  element  by 
the  X  •  1  element. 

the  negative  of  the  thrust,  such  as  from  a  propeller, 
that  the  vehicle  exerts  on  the  water. 


"x 


Vvf 

vvfp 

vhf 

vhfp 


the  velocity  of  the  c.g.  of  the  X  th  element. 

components  of  linear  and  angular  velocity  of  the 
vehicle  referred  to  axes  fixed  in  the  vehicle,  where 
i  -  1,2,...,  6. 

the  velocity  through  the  fluid  of  the  vertical  fins. 

the  projection  of  v^  on  the  |i5a  plane. 

the  velocity  threugh  the  fluid  of  the  horizontal  fins. 

the  projection  of  on  the  ^  §3  plane. 

weight  per  unit  length  of  the  X  th  element  in  a  vacuum. 

a  radius  vector  from  a  arbitrary  inertial  origin  to 
that  terminal  end  of  the  X  th  element  nearest  the 
towing  vehicle. 
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the  radius  vector  from  an  arbitrary  inertial  origin 
to  the  c.g.  of  the  vehicle. 

the  vector  from  the  towing  vehicle  to  the  towed 
vehicle. 

in  Appendix  A,  the  ’actor  from  the  c.g.  of  the  vehicle 
to  the  fin  axis. 

the  angular  deflection  of  the  vertical  fin  relative 
to  the  vehicle. 

the  angle  between  the  vertical  fin  and  the  projection 
on  the  5x5a  plane  of  its  velocity  through  the  water. 

the  angular  deflection  of  tH  horizontal  fin  relative 
to  the  vehicle. 

the  angle  between  the  horizontal  fin  and  the  projection 
on  the  plane  of  its  velocity  through  the  water. 

the  Kronecker  delta.  It  is  zero  unless  i  *  j  in  which 
case  it  is  one. 

coordinates  fixed  in  the  vehicle  with  their  origin 
at  its  c.g. 

the  index  denoting  the  cable  segment.  \  *  v  +  1 

denotes  the  vehicle, 

also  the  index  denoting  the  cable  segment. 


the  surface  integral  p  J  ^(dcp^/dnjdS,  cp^  is  the 
velocity  potential  due  to  motion  with  unit  velocity 
in  the  i  th  direction. 

the  number  of  segments  into  which  the  cable  is  divided, 
the  mass  density  of  the  fluid, 
the  angle  between  and  u^  . 

a  computational  abbreviation  for  +  FD^  +  FL^ 

+  (FIN)^,  where  (FIN)^  is  zero  unless  X  *  v  +  1. 


the  factor  relating  the  velocity  through  the  fluid 
of  the  X  th  element  to  the  normal  force  caused  by 
this  velocity. 

the  factor  relating  the  orientation  of  the  X  th 
element  and  the  normal  force  on  It  due  to  Its  velocity 
through  the  water. 

the  angular  velocity  of  the  X  th  element. 


STATEMENT  OF  THE  PROBLEM 


An  underwater  vehicle  is  towed  with  a  flexible  cable  by  a 
maneuvering  ship.  A  command  is  given  which  actuates  the  bow  planes 
of  the  vehicle  through  a  prescribed  angle.  The  vehicle  then  changes  its 
orientation  and  position.  It  is  required  to  determine  the  position  of 
the  towed  vehicle,  and  the  cable  tensions  as  time  varies. 

The  only  kinematic  boundary  condition  imposed  on  the  system  con¬ 
sisting  of  the  cable  and  underwater  vehicle  is  that  the  end  of  cable 
attached  to  the  towing  platform  follows  a  known  curvilinear  path  along 
which  the  components  of  velocity  and  acceleration  are  given.  It  is 
assumed  that  the  end  of  cable  attached  to  the  vehicle  if  extended,  would 
intersect  the  longitudinal  axis  of  the  vehicle.  This  occurs,  for  example, 
when  the  bail  is  an  arc  of  a  circle  whose  center  is  on  the  longitudinal 
axis. 


Other  items  which  are  taken  to  be  known  are  the  length,  mass  and 
diameter  of  the  circular  cable,  and  the  mass  and  dimensions  of  the 
vehicle. 


APPROXIMATE  PHYSICAL  SYSTEM 


The  cable  is  assumed  to  be  inextensible  (does  not  stretch)  but  its 
shape  is  flexible  in  form. 

The  length  of  the  cable  is  divided  into  v  integral  number  of  seg¬ 
ments  as  shown  in  Figure  1.  The  sections  are  imagined  to  be  short  rigid 
cylinders  with  universal  joints  at  the  junctions.  The  segments  are 
numbered  according  to  the  subscript  X.  X  =  1  for  the  segment  next  to  the 
towing  platform  and  X  **  v  +  1  is  the  underwater  vehicle.  The  assumption 
of  discrete  divisions  is  a  reasonable  one,  and  its  validity  improves 
as  v  increases.  The  difference  between  the  performance  of  the  actual 
and  assumed  cable  is  expected  to  be  small  as  v  becomes  large.  A  precedent 
for  this  assumption  appears  in  Reference  1. 
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In  the  problem  considered  the  cable  and  vehicle  do  not  roll.  This 
will  occur  if  either 

(1)  the  vehicle  is  roll  stabilized,  or 

(2)  the  shape  of  the  bail  on  the  vehicle  is  an  arc  of  a  circle  the 
center  of  which  is  on  the  longitudinal  axis  of  the  vehicle,  and  the 
cable  attachment  is  allowed  to  swivel  so  there  is  no  torque  transmitted 
from  the  laid  cable  to  the  vehicle. 

If  either  of  these  conditions  is  satisfied  then  roll  equal  zero  is 
a  solution  of  the  equations. 

The  vectors  in  the  following  problem  have  one  of  two  distinct 
notations:  (1)  the  customary  vector  bar,  or  (2)  a  subscript  indicating 
its  components.  The  first  method  will  be  used  when  possible. 

The  validity  of  a  vector  equation  is  independent  of  the  particular 
coordinates  in  which  it  is  expressed,  so  that  coordinates  are  generally 
left  unspecified.  If  the  vehicle  moves  through  an  undisturbed  fluid 
medium,  the  fixed  frame  may  be  taken  relative  to  the  earth.  However, 
if  the  fluid  medium  moves  with  a  velocity,  it  is  convenient  to  fix  the 
reference  frame  to  an  undisturbed  particle  moving  with  the  medium. 

Let  x^  be  the  radius  vector  from  the  origin  of  such  a  fixed  frame 

to  the  forward  terminal  point  of  the  \  th  element.  To  describe  the 
motion  of  the  system,  it  will  be  necessary  to  specify  the  motion  of 
each  cable  element  and  also  the  motion  of  the  underwater  vehicle.  The 
velocity  and  acceleration  of  each  of  the  terminal  points  relative  to 

JU 

the  inertial  frame  are  x  (t)  and  x (t).  The  position  of  the  bow  of  the 

\  K  _ 

underwater  vehicle  relative  to  the  ship  is  denoted  by  xv/g»  where 

x  =  x  . .  -  x,  .  It  is  not  easy  to  define  a  useful  quantity  "diversion" 
v/s  vrt 

while  the  towing  ship  is  maneuvering.  However  if  the  towing  ship  is  not 
maneuvering,  the  lateral  diversion  or  distance  can  be  shown  to  be: 


where  sq  is  a  unit  vector  pointing  from  the  bow  of  the  towing  ship 

toward  the  stern,  and  k  is  a  unit  vector  along  the  X3  axis  in  the 
inertial  frame,  as  shown  in  Figure  3.  The  rate  of  diversion  of  vehicle 
relative  to  ship  is  then 
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Further,  the  quantity  k  .  represents  the  position  of  bow  of  the 

vehicle  relative  to  the  sea  bottom,  providing  the_inejrtial  frame  is 
placed  on  the  sea  bottom.  Likewise  the  quantity  k  .  x^^  would  repre¬ 
sent  the  .ate  of  vertical  displacement  of  bow  of  vehicle  relative  to 
the  above-mentioned  frame.  The  tensions  in  t^e  cable  at  the  towing 

ship  and  at  the  towing  vehicle  are  T,  (t)  and  T  (t)  respectively,  and 

l  vri  _ 

a*-  th?  forward  end  of  X  th  element,  the  tension  is  (t).  In  a  similar 

manner  other  quantities  can  be  formulated.  These  quantities  are  ob¬ 
tained  from  a  solution  of  differential  equations  of  motion  of  the 
cable-vehicle  system.  The  formulation  of  the  equations  of  motion  is 
given  on  Page  23.  The  effect  of  the  accelerated  surrounding  fluid  on 
both  cable  and  vehicle  is  represented  more  accurately  in  the  con¬ 
siderations  originally  due  to  Kirchhoff  (Reference  3),  and  more  simply 
by  a  three-dimensional  generalization  of  the  expression  used  by  Walton 
and  Polachek  (Reference  1).  The  latter  method  is  employed  for  cable 
elements,  whereas  for  the  vehicle  two  approaches  are  considered.  The 
two  approaches,  designated  as  Approach  A  and  Approach  B,  are  described 
briefly  as  follows: 

1.  The  force  on  the  vehicle  is  calculated  on  the  basis  of  the 
same  assumptions  as  made  for  a  force  on  a  cable  element.  That  is, 

the  magnitude  is  proportional  to  the  volume  of  fluid  displaced,  and  if 
a  is  the  normal  component  of  the  linear  acceleration,  the  magnitude 
o?  the  force  is  proportional  to  |"an| ,  and  Its  direction  is  -'a^.  The 

proportionality  constant  for  the  vehicle  may  have  to  be  determined  ex¬ 
perimentally.  This  approach  has  the  undesirable  feature  that  trans¬ 
lation  and  rotation  are  not  coupled  through  the  reaction  of  the  ac¬ 
celerated  fluid. 

/ 

2.  The  force  on  the  vehicle  Is  given  by  hydrodynamic a 1  con¬ 
siderations  originally  due  to  Kirchhoff.  Physically  the  force  is  due 
to  the  accelerated  motion  generating  fluid  pressure  variations  over 
the  surface,  hence  translation  and  rotation  are  coupled  by  the 
motion  of  the  fluid.  This  approach  has  the  undesirable  feature  that 
it  complicates  a  formulation  already  lengthy. 

The  equations  of  motion  of  the  vehicle  using  either  Approach  A  or 
Approach  B  are  developed  separately.  The  resulting  systems  of  equations 
differ  from  each  other  considerably  in  form.  The  testing  of  the  two 
hypotheses  implies  a  large  experimental  program.  However,  both  are 
within  the  framework  of  a  reasonably  accurate  formulation  of  physical 
facts. 

To  determine  the  time  varying  quantities  such  as  diversion. 
Instantaneous  cable  tensions  and  path  of  vehicle,  we  begin  with  a  dis¬ 
cussion  of  the  dynamical  equations  of  motion  of  a  rigid  body.  The 
first  of  the  following  two  equations  describes  the  translation  of  mass 
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center  and  its  relationship  to  the  resultant  external  force.  The  second 
equation  describes  the  rotation  of  the  body,  and  specifically  it  relates 
the  resultant  external  moment  to  the  change  in  angular  momentum.  F  and 
M  with  suitable  subscripts  are  the  forces  and  moments  respectively  acting 
on  a  body,  H  is  the  angular  momentum  of  a  body  and  u  is  the  velocity  of 
the  mass  center  of  a  rigid  body.  The  two  equations  are: 

m^i/dt)  *  Fxl,  (\  *  1,2,...,  v  +  Ij  i  -  1,2,3) 
d^/dt  *  MXi,  (X  “  1,2,...,  v  +  1;  i  -  1,2,3) 


The  subscript  X  refers  to  the  X  th  segment  of  the  cable-body  system,  and 
the  index  i  refers  to  the  component  in  the  direction  of  the  axes  of  the 
fixed  frame  of  reference.  For  example,  F^  is  the  i  th  component  of  all 

forces  exerted  on  this  element.  The  reader  unaccustomed  to  the  use  of 
subscripts  to  denote  the  component  of  a  vector  should  refer  to  Appendix  B. 
While  both  equations  are  defined  relative  to  a  space  fixed  frame  x^ife, 
the  second  equation  can  also  be  written  with  respect  to  axes  fixed  to  the 
center  of  mass  of  any  body,  providing  all  quantities  entering  into  it  are 
referred  to  this  point.  The  force  F  may  be  expressed  a."  follows: 


Xi  AXi 


+  +  F,.. ,  +  (FIN) 


VXi 


GXi 


where  F.  is  the  sum  of  the  hydrodynamic  forces  which  depend  in  any  way 

on  the  accelerations  of  the  coordinates  of  the  end  points  of  the  elements. 
Py  ^  is,  with  the  exception  of  (FIN)^,  the  sum  of  all  hydrodynamic  forces 

which  depend  on  the  first  time  derivatives  but  not  on  the  second  time 
derivatives,  F  is  the  sum  of  the  forces  which  df-pend  only  on  the 

properties  of  the  element  and  on  its  orientation,  and  (FIN)^  is  the  force 

exerted  on  the  vehicle  due  to  the  deflection  of  the  fins  and  velocity 
of  the  water. 


THE  FORCE  F..  . 

AXi 

An  accelerated  vehicle  produces  accelerations  of  the  surrounding 
fluid  and  so  changes  its  momentum.  The  rate  of  change  of  linear  (or 
angular)  momentum  is  proportional  to  the  force  (or  moment)  producing 
this  change  which  in  turn  is  proportional  to  the  linear  (or  angular) 
acceleration  of  the  vehicle.  Ihus  the  vehicle  behaves  as  it  its  mass 
were  increased.  This  phenomenon  is  termed  "mass  accession"  or  added 
masses  and  added  mass  moments  of  inertia. 

In  Reference  1,  Walton  and  Polachek  assume  that  a  hydrodynamic 
force  is  proportional  to  that  component  of  acceleration  of  the  center 
of  gravity  of  the  cable  segment  which  is  normal  to  this  segment.  This 
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may  be  generalized  to  three  dimensions,  and  has  an  obvious  application 
to  a  long  thin  cylinder,  i.e.,  a  cable  segment.  However  when  we  come 
the  angular  momentum  equation,  there  is  no  allowance  made  for  the 
accompanying  added  mass  moment  of  inertia.  Approach  A  applies  this 
analysis  to  the  towed  vehicle  and  the  cable  segments. 

As  formulated  by  Kirchhoff  the  forces  are  linear  combinations  of 
the  accelerations.  The  coefficients  of  this  linear  form  are  surface 
integrals  of  the  velocity  potential  which  describe  the  fluid  flow.  How¬ 
ever,  some  of  the  coefficients  or  "aaded  masses"  may  be  obtained  experi¬ 
mentally.  Approach  B  applies  these  considerations  to  the  towed  vehicle. 


THE  FORCE  F^ 

There  are  two  hydrodynamic  forces  due  to  the  velocity  of  a  cable 
segment  through  the  water.  They  are  the  forces  normal  and  tangential 
to  the  cable.  each  r»f  theae  we  generalize  the  expressions  used  by 

Pode  (Reference  2)  to  three  dimensions.  The  assumptions  made  here  are 
as  follows: 

1.  The  tangential  force  has  the  direction  of  s  ,  where  s  is 

A  A 

discussed  in  Appendix  A. 


2.  The  magnitude  of  the  tangential  force  is  independent  of  the 
orientation  of  the  cable  segment. 

3.  The  magnitude  of  tangential  force  is  ( 1/2) . (constant) 

A  ■  Ci  uj  where  u^  is  the  velocity  of  the  segment  through  the  water, 
and  A  is  the  wetted  area  of  the  segment. 


4.  The  normal  force  is  in  the  plane  of  the  concurrent  vectors 

"a  and  u  ,  and  is  normal  to  "a  . 

AAA 

5.  l’he  magnitude  of  the  normal  force  is  the  product  of  two  factors 

X(Xl^)  where  cp^  is  the  angle  between  -"s^  and  u^.  Further 

t(<p^)  -  sin8  q>^,  and  X(^)  ■  (1/2)  p  ^  A  (constant)  -  QjU^.  (Refer  to 

Appendix  0  for  a  more  detailed  discussion.) 


THE  FORCE  F_, , 
v»Xl 

There  are  two  effects  that  contribute  to  F-, , .  The  first  is  the 

v»xi 

weight  of  the  cable  in  water.  The  second  is  composed  of  the  tensions 
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T^,  -  T\+i*  at  eac^  en<*  transmitted  by  adjacent  segments.  This  we 

suppose  to  be  propagated  instantaneously.  Since  action  is  equal  to 
reaction  and  since  segments  of  finite  length  will  have  finite  angles 
of  intersection,  this  force  will  not  be  along  the  axis  of  the  segment. 
A  diagram  of  these  forces  with  an  explanatory  notation  is  contained  in 
Appendix  C.  The  possibility  of  the  vehicle  exerting  thrust,  such  as 
from  a  propeller,  may  be  included  by  setting  -T^  *  -Tp8^j,  where  Tp 

is  the  magnitude  of  the  propeller  thrust. 


THE  FORCE  (FIN) t 


The  lift  and  "drag"  components  of  the  resultant  hydrodynamic  force 
on  the  fins  are  assumed  to  follow  quasi-steady  and  real-flow  formula¬ 
tions.  Forces  which  are  zero  when  the  rate  or  rotation  of  the  fins  is 
zero  are  neglected.  Th*  fins  on  the  vehicle  are  in  both  the  vertical 
and  horizontal  planes. 


Let  v  ^  denote  the  velocity  of  the  vertical  fins  through  the  water, 
and  v^  denote  the  projection  of  v^  on  the  §i5s  plane.  (Refer  to 
Appendix  F).  The  following  assumptions  are  made: 


1.  The  force  is  called  lift  if  it  is  normal  to  the  plane  con¬ 
taining  both  vvfp  and  n;  and  a  force  is  called  drag  if  it  is  in  the 

direction  -v  ,  .  Lift  is  positive  when  it  tends  to  increase  the  x  - 

vxp 

coordinate  from  its  initial  steady  state  value. 

2.  Both  of  the  forces  are  proportional  to  lvv£p|2* 


3.  Each  force  is  a  function  of  the  angle  of  attack;  where  the 
lift  is  an  odd  function  and  the  drag  an  even  function  of  the  angle  of 
attack. 


4.  The  factors  of  proportionality  in  these  forces  are  r3  and  C4 
for  lift  and  drag,  respectively.  Then  C3  and  C4  include  variations 
due  to  fin  area  and  variations  in  and  cDs  which  are  coefficients  of 
lift  and  drag  respectively. 

The  forces  on  the  vehicle  due  to  horizontal  fins  are  considered  in 
exactly  the  same  way.  The  lift  due  to  these  fins  is  normal  to  v^ 

(the  velocity  through  the  water  of  the  horizontal  fins) ,  and  is  it*  the 
vertical  plane.  The  horizontal  and  vertical  fins  are  taken  to  be  at  the 
same  distance  from  the  bow,  although  this  assumption  is  not  necessary 
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and  may  be  dropped  for  different  vehicle  designs.  However  this 
assumption  implies  v^^  *  vyj  .  For  this  pair  of  fins  we  obtain  two 

other  constants,  C5  and  Cg ,  for  the  lift  and  drag,  respectively.  Also 
the  lift  is  considered  to  be  positive  when  it  tends  to  increase  the 
x „  coordinate  from  its  initial  steady  state  value.  A  detailed 

VT 1  i  j 

treatment  is  found  in  Appendix  F. 

Now  consider  the  second  set  of  equations  of  motion  which  relate 
the  changes  of  angular  moments  of  the  vehicle  or  of  the  cable  elements 
to  the  external  applied  moments.  These  equations  are: 

«u-V 


In  all  cases  the  angular  momentum  and  moment  of  the  forces  are  taken 
about  the  center  of  gravity  of  the  segment  involved.  As  shown  in 
Appendix  B, 


u  - 

"Xi 


„  -X  *  T\ 
"Xi"  jj  '  Vji 


e.  u),  ,u),  I" 
1st  Xj  Xs  jt 


where  uj 


Xj 


is  the  angular  velocity  of  the  cable  element  and  I 


X 

ij 


The  quantity  is  the  i  th  component  of  the  vector  from  the  c.g.  to  a 

typical  small  mass  m  sud  the  SwwwCiuu  Is  carried  out  over  the  X  th 
element.  The  above  equation  is  referred  to  a  frame  whose  origin  is 
attached  to  the  center  of  gravity  and  whose  axes  are  parallel  to  the 
axes  of  the  Inertial  frame. 


To  apply  the  above  equation  to  the  vehicle,  it  is  more  convenient 
to  refer  the  rate  of  change  of  moment  of  momentum  to  a  set  of  axes 
imbedded  in  the  vehicle.  Let  the  superscript  prime  denote  the  quantities 
in  this  frame,  then  the  above  relation  becomes: 


H' . ,  . 
v+l,i 


v+l.t  JJ 


r.v+1  _ 


.  |V+1 


V*1.J  Ji 


>  (I)'  (l)1 

"ist  v+l,j  v+l,s 


T.V+1 

jt 


To  obtain  o',  I'  and  the  H'  relative  to  the  body  fixed  frame,  let  the 
rotation  of  the  vehicle  relative  to  the  fixed  inertial  frame  be 
described  by  a  transformation  of  coordinates  as  follows: 


cij5j  °r  5k 


cikxi 


The  c^j  are  obtained  in  terms  of  previously  defined  quantities  as 
shown  in  Appendix  G.  Since 
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h  . ,  .  *  c .  a' 

v+l,i  in  v+l,n 


and 


"vtl,i  cinVn,n’ 


“*  v+l,j  “  Vl,AC/j 


0)'  ,  -  0)  ,  ,  .C..,  to’  ■  U)  ,  c, 

v+l,j  v+1,4  ty  vfl,s  vfl.k  ks 

the  preceding  expression  for  H  ,  relative  to  the  body  fixed  frame  now 

VtI  1 1 

takes  the  following  form: 

•  •  tv+l  *  . v+1 

Hv+l,i  "  u)v+l,iIjj  “  u)v+l,icinCijIjn 


„  „  „  T.vfl 

nst  v+l,£wv+l,k“in”£j”ks“jt 


-  r>  n\ 


The  products  and  moments  of  inertia  relative  to  the  body  fixed  frame 
are  treated  in  Appendix  B. 


The 


value  of  for  X  ■  l,2,...,v  (i.e.,  for  cable  segments)  is 


more  readily  obtained,  since  the  1^  depend  on  by  a  known  functional 

relation.  In  fact  they  finally  simplify  to: 

♦  • 

_  /  i2  / 1  f\\  ... 


TY  -  /  **  *  Y  ... 

l\t  ~  '“W 


Xi 


The  right  side  of  the  angular  momentum  equations  contains  the 
applied  external  moments  with  respect  to  the  c.g.  of  the  body.  For 
convenience  we  express  these  in  a  manner  similar  to  that  in  which  the 
forces  were  decomposed. 


"x  ■  “gxi  +  V  +  mau  + 


where: 


is  the  moment  of  those  forces  which  depend  on  neither 
the  velocities  nor  accelerations, 

My  .  is  the  hydrodynamic  moment  of  those  forces  other  than 
*  those  due  to  the  fins  depending  on  the  velocities  but 
not  on  the  acceleration, 


M 


AXi 


is  the  hydrodynamic  moment  of  all  forces  dependent  on 
acceleration, 


(FINM)^  occurs  only  for  the  vehicle  and  is  the  moment  due  to  the  fins. 
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THE  MOMENT  M. 


GU 


For  a  cable  element  this  moment  is  due  only  to  forces  exerted  by 
adjacent  elements.  For  the  vehicle  we  recall  that  in  some  designs 
there  may  be  a  propoller  thrust  whose  line  of  action  does  not  pass 
through  the  c.g.  of  the  vehicle.  Should  this  be  the  case,  the  ad¬ 
ditional  moment  of  the  propeller  thrust  about  the  c.g.  must  be  included. 
Furthermore  there  occurs  an  external  moment  on  the  vehicle  whenever  the 
center  of  gravity  and  center  of  buoyancy  do  not  coincide.  Appendix  C 
may  be  consulted  for  details. 


THE  MOMENT  M^ 

The  physical  assumptions  relating  the  forces  due  to  the  velocity 
have  already  been  listed.  Since  translation  and  rotation  of  a  cable 
element  occur  simultaneously,  one  end  travels  through  the  water  faster 
than  the  other.  This  motion  creates  an  opposing  hydrodynamic  moment. 
The  computation  of  this  moment  is  carried  out  in  Appendix  D.  The 
moment  of  forces  due  to  velocity  on  the  vehicle  is  treated  in  an  en¬ 
tirely  similar  manner. 


THE  MOMENT  MAU 

The  hydrodynamic  moment  is  zero  according  to  the  hypothesis  stated 
in  Reference  1.  This  holds  for  all  elements  designated  by  X  ■  l,l,...,v, 
and  also  for  the  vehicle  designated  by  X  ■  v  +  1  in  Approach  A.  However 
for  the  formulation  discussed  in  Approach  B,  there  is  a  moment  on  the 
vehicle  due  to  the  acceleration  of  the  medium.  In  this  instance  tin 
moment  is  calculated  in  exactly  the  same  manner  as  the  previously  dis¬ 
cussed  hydrodynamic  force.  The  detai  are  found  in  Appendix  E. 


THE  MOMENT  (FINM>1 

The  moment  acts  only  on  the  vehicle.  The  fins  are  symmetric  airfoils. 
(FINM^  is  the  moment  about  the  c.g.  due  to  the  lift  and  drag  on  the  two 

pairs  of  fins.  The  exact  expression  is  contained  in  Appendix  F. 

We  now  consider  the  mathematical  approach  to  the  problem. 
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THE  MATHEMATICAL  PROBLEM 


The  equations  to  be  solved  are: 

m^du^dt)  -  Fu,  (X  -  1,2 . v  +  1;  i  -  1,2,3) 

dl^/dt  -  Mxi,  (X  -  1,2 . v  +  1;  i  -  1,2,3) 

where  du^/dt  is  the  linear  acceleration  of  c.m.  of  the  X  th  element, 

and  dH^/dt  is  an  expression  linear  in  the  angular  accelerations  of  the 

X  th  element.  In  connection  with  these  equations  there  ere  two  funda¬ 
mental  considerations. 


1.  The  right  side  of  »ny  equation  of  the  first  set  contains 

T  -  T  +  F*  ,  and  the  right  side  of  any  equation  of  the  second 
X  *  X  '  XI 

set,  except  for  X  *  v  +  1,  may  be  written  “^j^8  ij^Xk  +  TX+l,k^  +  “xV 
where  F*  and  M*  are  expressions  which  do  not  contain  the  terms  T  .  By 

summing  the  first  set  of  equations  from  X  ■■  v  +  1  to  X  *  n,  an  expression 

for  T  .  is  obtained.  Thus  T  .  +  T  .  .  (or  T  +  T,  , ,  , )  may  be  obtained. 

p.t  M-i  M-+l,i  Xi  X+l,i  J 

When  this  is  substituted  in  the  equation  for  angular  momentum  all  T^  are 

eliminated,  and  a  set  of  equations  in  fc  and  'k  is  obtained.  This  set  can 

Ko  int-oorat-pd  using  numerical  orocesses. 


2.  The  second  fundamental  consideration  is  that  there  are  constraints 
among  the  elements.  Suppose  we  make  the  natural  choice  and  take  as 
coordinates  the  position  of  the  ends  of  the  elements.  We  automatically 
have  a  connected  chain  of  elements.  The  constraints  are  that  each  element 
has  a  given  constant  length. 


The  position  of  a  rigid  body  is  completely  determined  by  the  location 
of  any  three  of  its  points,  provided  that  these  points  are  noncolinear. 

The  distances  between  these  points  do  not  change.  Hence  there  are  three 
constraints  among  the  nine  coordinates  and  six  degrees  of  freedom. 


In  the  problem  under  consideration  suppose  for  the  moment  that  the 
ends  of  the  rigid  element  are  not  connected  to  the  adjacent  elements. 
Suppose  the  coordinates  of  the  two  end  points  of  the  element  are  specified. 
Because  there  is  one  constraint  among  the  six  coordinates,  five  degrees 
of  freedom  remain.  When  the  location  of  either  vehicle  or  cable  element 
is  specified  by  giving  thr.  positions  of  its  end  points,  the  degree  of 
freedom  of  rotation  about  a  line  joining  the  two  points  is  lost  because 


roll  is  assumed  to  be  zero.  The  condition  of  zero  roll,  id. 


Xj  ‘  Xj 


=  0, 
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yields  a  relation  among  the  three  angular  velocities  referred  to  the 
initial  system.  Further,  there  will  be  a  relation  among  the  three  com¬ 
ponents  of  angular  acceleration.  Then  for  any  element  any  component  of 
the  angular  velocity  may  be  found  in  terms  of  the  other  two.  This  is 
also  true  for  the  angular  acceleration,  and  therefore  also  true  ::or 
the  angular  momentum  and  its  time  derivative. 

When,  as  in  the  preceding,  the  roll  about  the  longitudinal  axis 
is  zero,  the  angular  velocity  may  be  readily  determined  in  terms  of  the 
positions  of  the  end  points  of  the  element  and  their  derivatives.  The 
expression  obtained  for  u>^  is  dependent  upon  the  fact  that  the  element 

has  a  constant  length  (c.f.  Appendix  A);  i.e.,  d(length  of  element)/dt  -  0. 
Similarly  the  relationship  of  <1)^  to  the  coordinates  of  the  end  points  and 

their  derivatives  is  dependent  upon  the  constant  length  of  the  element. 
Therefore  there  is  a  dependency  relation  amo.jg  the  three  equations  of 
angular  momentum  and  the  condition  that  the  length  of  the  element  is 
constant. 

The  preceding  discussion  was  under  the  hypothetical  condition  that 
the  ends  of  the  element  were  free,  which  implied  five  degrees  of  freedom. 

In  our  problem  the  ends  of  the  element  are  not  free.  The  forward  end 
of  the  X  +  1  st  element  must  coincide  with  the  after  end  of  the  X  th 
element.  This  eliminates  three  degrees  of  freedom  for  the  X  +  1  st 
element,  so  that  each  element  has  two  degrees  of  freedom. 

Fcr  example  the  element  next  to  the  towing  vehicle  has  a  given 
^  t-ha  nnsIHmt  of  If*  forward  .-nd  <r  Rnacified.  Clearly  the 

position  of  its  after  end  lies  on  a  sphere  centered  about  the  forward 
end,  and  may  be  specified,  for  Instance,  by  giving  its  latitude  and 
longitude.  When  these  two  values  are  specified  for  the  first  element, 
two  more  must  be  given  to  determine  the  position  of  the  second  element, 
and  so  forth. 

Thus  there  are  two  degrees  of  freedom  for  each  element  and  three 
angular  momentum  equations,  of  which  only  two  are  independent  when  account 
is  taken  of  the  constancy  of  length  of  an  element. 

We  may  use  these  two  angular  momentum  equations  for  the  two  inde¬ 
pendent  angular  velocities  and  obtain  the  third  angular  velocity  from 
the  constraint.  On  the  other  hand,  anv  nonindependent  equation  may  be 
dropped  from  a  dependent  system.  To  drop  the  constraint  it  is  sufficient 
to  know  that  the  three  remaining  equations  are  not  independent  of  the 
constraint.  As  mentioned  previously  and  as  shown  in  Appendix  A,  the  ex¬ 
pressions  obtained  for  and  <1^  are  dependent  on  the  constraint,  so 

that  instead  of  using  two  angular  momentum  equations  and  the  constraint, 
we  may  use  the  three  angular  momentum  equations. 
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DISCUSSION  OF  METHOD  OF  SOLUTION  OF  EQUATION 
OF  MOTION 


Having  taken  as  coordinates  the  end  points  of  the  elements,  the 
cu's  and  u's  which  occur  in  the  equations  of  motion  are  obtainable  in 
terms  of  derivatives  of  these  coordinates.  The  3(vH)  linear  momentum 
equations  determine  the  3(v+l)  forces  T  ,  and  the  3(v+l)  angular  momen¬ 
tum  equations  then  determine  the  3(v+l)A  coordinates  of  the  end  points 
of  the  v  +  1  elements. 

The  final  3(v+l)  equations  turn  out  to  be  linear  in  the  linear 
accelerations.  At  this  point  we  anticipate  using  a  subroutine  for 
solving  differential  equations  such  as  SHARE  program  GL  AIDE  1  (Ref¬ 
erence  4).  This  necessitates  a  sequence  of  availability  of  information 
which  is  obtained  by  using  previous  values  of  positions  or  velocities  to 
obtain  the  present  value  of  acceleration.  Thus  we  use  Known  values  of 
positions  and  velocities  to  compute  the  coefficients  of  the  linear 
equations  for  the  accelerations  and  the  constants  on  the  right  sides  of 
these  equations.  The  linear  system  may  then  be  solved  numerically  to 
obtain  values  for  the  acceleration.  This  satisfies  the  requirements  of 
the  numerical  program. 


EQUATIONS  OF  MOTION  FOR  APPROACH  A  AND  APPROACH  B 


In  this  part  the  expressions  for  forces  and  moments  which  act  on 
the  vehicle  are  substituted  into  the  equations  of  motion.  The  final 
differential  equations  of  motion  for  Approach  A  are  denoted 
Equations  (1)  and  (2),  and  those  for  Approach  B  are  denoted  by 
Equations  (3)  and  (4) .  A  solution  of  either  set  of  equations  would 
yield  the  components  of  velocity  of  the  displacement  coordinates  at  each 
of  the  cable  joints  and  also  at  the  bow  of  vehicle. 

Both  formulations  are  now  presented.  The  first  formulation  con¬ 
siders  the  assumptions  of  Approach  A,  and  the  second  formulation 
considers  the  assumptions  of  Approach  B. 
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The  Equitlooe  of  Motion  for  *«roech  l 


31 


UttO 


32 


for  the  vehicle: 


•11  UCM  containing  accelerations  are  written  cm  the  l"ft  side,  and  all  other  tens  aie  wrltta-  on  the  right  side.  The  result  is: 


proceeding  syateia  contains  3v  +  3  equations.  The  3(w4-l)  unknown  are  x-  ,  x  ...  x  -  .  The  system  is  linear  in  the  accelerations  and 


CONCLUSIONS 


The  equations  of  motion  for  a  towed  underwater  vehicle  have  been 
derived.  The  cable  was  imagined  to  consist  of  a  finite  number  of  rigid 
cylinders.  Two  hypotheses  were  made  relating  the  hydrodynamic  force 
and  moment  on  the  vehicle  and  its  acceleration,  and  the  equations  were 
developed  separately. 

The  equations  of  motion  are  involved,  the  most  desirable  of  the 
two  hypotheses  above  being  the  most  intricate.  In  these  equations 
there  occurs  a  sevenfold  sum.  Complexity  of  the  equations  is  the 
factor  which  will  limit  the  computability. 

The  principal  restriction  in  the  equations  of  motion  is  that  the 
vehicle  does  not  roll,  nor  does  the  cable.  On  the  other  hand,  the 
equations  are  general  in  that  both  the  towing  and  towed  vehicle  may 
move  in  a  highly  arbitrary  fashion. 
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(Reverae  Page  36  Blank) 


APPENDIX  A 


GEOMETRY 


This  appendix  contains,  with  the  exception  of  the  coordinate 
transformation  in  Appendix  G,  all  of  the  diverse  kinematical  considera¬ 
tions  relevant  to  the  problem. 


0  (Origin  of  Inertial  Frame) 

FIGURE  Al.  THE  X  th  CABLE  SEGMENT 

Let  s ^  be  a  unit  vector  tangent  to  the  cable  segment  of  length  d^. 
Referring  to  Figure  Al,  and  with  the  rules  of  vector  addition,  we  write 

Vi  ■  \  +  \V  °r  \  ■  <\+i  -  v/dx 

Since  the  length  of  cable  is  constant,  then  d(l/d^)/dt  *  0,  and  the 
derivative  of  s^  with  respect  to  time  gives: 

\  ■  <Vi  ‘  *i>dx  <A'1> 

In  i  similar  manner: 

“x  •  <Vi  *  V/dx 


-37- 


Since  the  center  of  mass  and  center  of  volume  of  a  cable  segment 
coircide,  by  the  rules  of  kinematics,  we  obtain 

■A.  mmm  __ 

%  ’  xc» +  v<V2)\ 

Vl  ■  xc.  +  “x  x 

Subtracting  these  two  equations  we  get 

\  -  Vi  ■  “x  x  Vx 

where  is  the  angular  velocity  of  the  cable  segment.  The  vector  product 

with  respect  to  a  gives: 

A 

•x  x  (xx+i  •  V  ■  Vx  x("  x  V  ■  Vx'V5-V*x 

Since  the  cable  segment  is  assumed  not  to  roll,  then  u)  .8  >0  . 

_  A  A 

Solving  the  previous  equation  for  tu ,  and  with  the  aid  of  Equation  A-l, 

—  —  A.  A.  JLA 

we  get  u).  ■  8  x  s  .  Since  s  x  s  “  0,  then  it  follows  that 

a  A  u  \  A  AA 

“x  M  sx  x  SX  • 

By  referring  to  Figure  A2,  the  location  of  the  center  of  mass  of 
vehicle,  with  respect  to  the  inertial  frame  may  be  expressed  in  the 
1  ** *> 

*w*.*wn*  w  • 

XG  “  (dv+l  "  dv+2^xv+l^dv+l  +(dv+2/dv+l)xvf2 

Here  the  velocity  and  acceleration  of_the  center  of  mass  are  the  usual 
first  and  second  time  derivatives  of  x  , 

VJ 

The  important  unit  vector  q  ,  as  shown  in  Figure  A3,  has  the 
following  properties:  (a)  it  is  in  the  plane  of  u^  and  s  ,  (b)  it  is 
normal  to  s^  ,  (c)  it  is  in  the  sense  that  the  direction  between 
q^  and  -u^  is  less  than  n/2,  and  (d)  it  has  unit  length.  The  analytical 
statement  of  these  conditions  is: 

1.  q^  ■  ou^  +  @s^,  for  some  constants  a  and  f}  . 

2-  \'a\  *  0  »  ('3)  <  °»  *  1 
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FIGURE  A2.  MEASUREMENTS  ON  VERTICAL  PROJECTION  OF  VEHICLE 


“x 

FIGURE  A3'.  UNIT  VECTORS  q>  ,  s  ,  AND  VELOCITY 
OF  CABLE  SEGMENT 

When  one  takes  the  dot  product  of  both  sides  of  (i)  with  s  ,  the  left 

A 

side  is  zero  by  (2) ;  then 

au  .  s  +  P  =  0 
A  A 
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If  this  value  of  0  is  substituted  in  (1),  the  result  is 

\  ’  “K  '  Y  Vx1  '  “K. 11  <\  *  Tx)] 

From  condition  (4)  and  with  this  expression  for  q  ,  we  obtain 

A 

1  ■  “x  x  <“x  x  Y-X  *  (\  *  V 

■  OpS^  .  <U^  X  s^)  X  [s^  X  (U^  X  8X)] 

■  \  ■  [(“x  -  V  •  V 

J- 

.  s  *  1,  and  solving  for  a,  where  y  is  either  ±  1,  we  obtain 

A. 

a  ■  y/lu^  x  8XI 

To  determine  whether  y  -  +  1  or  y  -  -1  is  applicable,  condition  (3)  may 
be  used.  Applying  this  condition,  the  following  result  is  obtained: 

\  ■  \  'Y  “x  •  \  *  <“x  x  \)l  «  «xl  <0 

Since  the  quantity  after  the  symbol  y  is  positive,  it  follows  that 
y  ■  -1.  Therefore 

\  *  *‘x  11  *  V'K  *  »xl 

Vehicle  Geometry 

d  , ,  is  the  distance  from  the  cable  attachment  to  the  tail, 

d^2  is  the  distance  from  the  cable  attachment  to  the  c.g., 

d  is  the  distance  from  the  cable  attachment  to  the 

v+3 

effective  axis  of  the  horizontal  fins. 

d  is  the  distance  from  the  cable  attachment  to  the 
v+4 

effective  axis  of  the  vertical  fins, 

d  , .  is  the  distance  from  the  cable  attachment  to  the 
v+5 

nose  (dv+5  £  0), 

d  is  the  distance  from  the  cable  attachment  to  the 
v+6 

center  of  buoyancy. 


Since  s 

A 
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Let  x  be  the  vector  from  the  c.g.  to  the  fin  axis 


has  the  sense  of  a  vector  from  the  cable  attachment  to  the 
c.g.  If  the  d's  are  positive  measured  in  the  sense  of  s^,  for  any 
placement  of  the  c.g.,  fins,  and  cable  attachment,  then: 

x  *  (d^+3  -  d^)  s^,  (horizontal  fins) 

*  *  <dv+4  "  dv+2}  W  (vertical  flns> 

Summary  of  results 

1.  The  linear  velocity  and  acceleratin  of  the  center  of  \  th 
segment  of  cable  with  respect  to  tne  inertial  frame  are: 

(vi +  v/2  ■  <vi +  \)n 

2.  The  linear  velocity  and  acceleration  of  the  vehicle  with 
respect  to  the  inertial  frame  are: 

VflVl  +  b  Xv+2 

XG  "  a  Xv+1  +  b  Xv+2 

where  a  -  (d^  -  d^Vl  ’  b  “  dNH-2/dv+l 

3.  The  vectors  s  ,  <#  ,  and  x's  are  relat  follows: 

A  A 

*x  ■  <Vi  •  VV  \  ‘  <VrVV  *x  ■  <Vi  •  Vdx 


\  •  'Xx  x  <“x  x  V7|"x  x  “x1 


(Reverse  Page  42  Blank) 
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In  this  appendix,  the  linear  and  angular  moments  and  their 
derivatives  of  the  X  th  element  are  considered. 


Consider  a  small  element  of  mass  ny  as  shown  in  Figure  Bl.  Since 
originates  fro r  the  c.m.  of  the  cable  element,  it  follows  from  the 
definition  of  center  of  mass  of  bodies  that 

1  =  °»  and  I  “c^i  =  ° 


Let  is.  be  the  components  of  angular  velocity  of  the  \  th  cable 

A 1 

element;  then  the  velocity  and  acceleration  of  a  typical  small  element 
of  mass  are: 


vel  of  m^  = 


111 


+  e  .  . .  a),  q 

iJK  Xj 


acc  of  m 


(Vl.i  +  ^ 


+  eij<^a>XjPCk  +  eijka,\  jPCk 


In  the  above  expressions,  the  summation  convention  tias  been  used. 
Within  one  term,  no  small  latin  index  is  to  be  repeated  more  than  once. 
When  it  is  repeated  once  within  a  term  a  summation  of  this  index  from 
1  to  3  is  implicitly  understood.  If  a  small  latin  index  is  not  repeated 
the  equation  is  to  hold  for  that  index  having  each  of  the  values  1,  2,  3. 
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The  convention  does  not  apply  to  Greek  indices,  and  may  be  suspended 
by  explicitly  stating  so.  Further,  eijjc  is  the  conventional  three- 

index  symbol,  where  i  ■  1,2,  or  3,  and  j  and  k  have  the  same  range.  If 
any  two  indices  are  equal,  then  e  ■  0.  If  i,  j,  k  are  an  even 

permutation  of  1,  2,  3  the  symbol  has  the  value  of  +  1.  If  i,  j,  k  are 
an  odd  permutation,  its  value  is  -1.  For  example,  the  i  th  component 
of  I  x  b  le 

The  rate  of  change  of  linear  momentum  of  the  X  th  cable  element  is: 
d£(Linear  momentum)^^  /dt  ■  d£  ^  ny  (velocityK^J  /dt 


+  ®ijk  “\j  L  mCpCk  +  eijk(UXj  I  mCPCk 


where  the  sums  are  carried  out  as  £  ranges  over  the  X  th  element.  The 
above  expression  reduces  to  the  following: 


dj^  (Linear  momentum)^ jJ  /dt  ■  (mass^) 


+ 


Lee  be  the  anguiar  momentum  about  the  c.g.,  then  the  rate  ot  change 
of  the  angular  momentum  of  the  \  th  cable  element  is: 


d<Hu)/dt  -  l  01Jk  PCJ  <„cc)Ck  -  ixl 


"u  '  eiJk  i+1'2  ’Sk  l  "CPCJ  +  ei)k  I  VtApq1 


U).  pr 

i  xp  Cq 


+  *tjk  L  ekpq“xp  (eqPt“xSkCt)  > 

in  which  «  ^st^s^k’  and  ^acc^k  derived  earlier  in  this  section 

have  been  substituted.  For  further  simplification  we  recall  the 
following  relation  (Reference  5): 

e,  ,.e.  »  e,.,e,  *  6*  6d  -  6*  6^ 

ijk  kpq  kij  kpq  p  q  q  p 
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so  that  the  second  term  in  the  above  expression  for  may  be  written 
as  follows: 


euk  I  VuViVc-i '  <6p  6q  ‘  6i  6P  I  Wc«% 

‘  V  I  mCpCjpCJ  '  "kj  I  "CPC1PC3 
■  “u  *jj  '  »Xj  Ilj  ’ 

wh««  I»j  -  I  »cPCjP(;j  «nd  IXM  -  I  »cpcl„a 


The  contribution  of  the  last  term  in  may  be  rearranged  into  the 

following  form,  where  £  ranges  over  the  X  tb  element  and  is  indicated 
by  the  summation  symbol: 


eijkekpqCqst 


Z  Vli’ltVu 


9 


which  reduces  to  the  following  form: 

(51  6^  -  6*  6^)  e  ,u>,  u),  Y  nup.  .p,.^ 

p  q  q  p  qst  xp  \s  L  i;HCJpCt 

=  <ti  m  <»  -  e  (l)  (ix 

Xi  kb  jst  jt  1st  X J  \a  jt 

X 

jt 

in  j  and  t  implies  that  the  first  term  of  this  sum  is  zero,  which  is 
proved  as  follows: 


where  as  before  I 


jt 


z 


"CpCJpCt 


Moreover  the  symmetry  of  I 


2e,  AA  *  e.  tu, .1*  +  e,  „u>,  .u>  I* 
jst  Xi  Xs  jt  jst  \t  jt  jst  Xi  Xs  tj 


<D.  . 0).  I*  (e.  -  e.  .) 

Xi  X8  jt  jst  jst' 


Finally  the  expression  for  for  the  X  th  cable  element  becomes: 

H.  .  *  w,  ,1*  "  ,lY  •  e,  .m ,  .u>  I* 

Xi  Xi  jj  Xj  ji  ist  Xj  Xs  jt 

X 

where  Ijt  is  summed  over  the  X  th  element.  This  expression  for 
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holds  whenever  the  expression  used  for  (acc)^  gives  the  acceleration 

of  the  C  th  particle  element  in  the  inertial  frame.  Since  the  latter 
is  a  vector,  it  is  valid  both  in  the  inertial  and  in  a  frame  fixed  in 
the  body.  However,  enters  into  the  equations  of  motion  in  moving 

coordinates  directly  only  when  the  origin  of  the  moving  frame  is  at  c.g. 

To  apply  the  H  expression  to  the  vehicle,  it  is  desirable 
v+i 

to  express  1^  as  a  time  independent  quantity.  With  this  in  mind,  we 

imagine  a  frame  fixed  in  the  vehicle  at  its  c.g.,  consequently  the 
expression  ^  with  reference  to  this  frame  is 


H' 


uv 


v+l,n  v+l,n  jj 
However 


i!*1  -  «. 


N*1,J 


IlVH-e  -• 


.,v+l 


jn 


nst  vH,j  v+1,8  jt 


H'  ,  , 

v+l,n 


0)' 


v+l,n 
c.  _H' 


v+l,n 


are  vectors,  so  that 


Hv+l,i  *  cinHv+l,n  *  %-l,i  “  cin“vfl,n 


®vfl,J  *  ®vfl,AciJ*  ®vfrl,J  "  ®vU,XC4j’  ®v+l,s  "  ®v+l,kck3 
where  the  c^  are  the  coefficients  found  in  Appendix  G, 


Finally  a  substitution  of  the  above  quantities  into  the  expression  for 
HvH,i  glve8! 


H 


U) 


vH-1,1  V41,i*jj 


.I'.Y*1  -  «  . .  ,c.  C„  .1 


t  v+1 
\H-l,i''invijJ’jn 


i  v+1 

”  enstU)v+l,£U,\H-l,kcinC-ejCksIjt 
•  • 

where  the  H,  the  u>'s,  and  u>'s  are  rt  lative  to  the  inertial  frame  and  the 
I's  are  calculated  relative  to  the  vehicle  fixed  frame  of  reference. 

The  expression  for  a  cable  element,  however,  may  be  expressed 

with  reference  to  the  inertial  frame.  The  moments  of  inertia  with 
respect  to  the  inertial  frame  ares  (where  A  *  nyd3 / 12) 
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2Q  r*se*rfm±  tzv  v  v  05* 


As  s  , 

y  >’ 


As  s  , 

y  ** 


As  s 
z  z 


where  sx>  s^,  sz  are  components  of  the  unit  vector  tangent  to  the  cable 

vector,  and  these  components  are  known  at  any  instant  of  time.  Here  the 
cable  diameter  d  is  assumed  to  be  made  much  less  than  the  length  of  a 


cable  element. 


APPENDIX  C 


GEOMETRICAL  FORCES 


There  are  two  types  of  geometrical  forces  on  each  cable  element. 
They  are  the  weight  of  cable  in  water,  and  the  force  exerted  by  ad¬ 
joining  elements.  The  vehicle  possesses  these  forces  also,  and  in 
addition  has  a  moment  about  the  c.g.  due  to  the  displacement  of  the 
center  of  buoyancy  from  the  c.g. 


\i 


-X. 


X  + 


T 


Xi 


FIGURE  Cl.  CABLE  TENSIONS  AND  NET  WEIGHT 


Let  s^  be  the  unit  tangent  vector,  w^  the  weight  per  unit  length  of 
cable  in  air,  and  b^  the  buoyant  force  per  unit  length  of  cable. 

The  force  on  a  cable  element  is  shown  in  Figure  Cl: 


FGXi  =  TXi  *  TX+l,i  +  (WX  "  \)dX631 
and  the  moment  of  these  forces  about  the  c.m.  of  the  element  is: 


MGXt  “  (dX/2)  l'"eijksXjTXk  +  eijksXj("  TX+l,k^ 


"  '  *  eijkdX8Xj(TXk  '  Vi,^ 

The  force  and  moment  on  the  vehicle  are: 

FGvfl,i  =  Tv+l,i  ’  Tv+2,i  +  <wv+l  '  bv4-i)dv+l631 
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MGv+l,i“  eijksv+l,j^dv+2  ‘  dv+l^  Tv+2,k  "  dv+2eijksv+l, jT\H-l,k 
"  eijksv+l,j63k(dv+6  '  dv+2)bv+ldv+l 

where  =  (buoyant  force  on  vehicle) /d^^.  Note  that  the  force 
applied  to  the  stern  of  the  vehicle  is  (-T^^  jK 
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APPENDIX  D 


VELOCITY  FORCES 

There  are  two  forces  dependent  on  the  velocity}  the  tangential 
and  normal  forces. 

We  assume  as  in  Reference  2  that 

1.  The  tangential  force  has  the  directions  s^. 

2.  The  magnitude  of  the  tangential  force  does  not  depend 
on  the  orientation  of  element. 

3.  The  magnitude  of  the  tangential  force  is 
(l/2)pt^A(constant)  ■  C^d^,  where  u^  is  the  velocity 
through  the  water. 

4.  The  normal  force  is  in  the  direction  q^,  that  is  in  the 
plane  of  the  concurrent  vectors  1s^  and  u^  and  normal 

t0  V 

5.  The  magnitude  of  the  normal  force  is  the  product  of  two 

factors  where  ^  *s  the  angle  between  -s^  and  u^. 

In  addition 

U\)  *  8  inS  \  * 

\(^)  -  (l/2)pvl^A(constant)  *  C^d^ 

The  u^  appearing  above  varies  along  the  element  due  to  its 
rotation.  We  ask  what  error  is  committed  using  the  u^  of  the  midpoint 
of  the  cable  element. 

Let  A F  be  the  magnitude  of  some  force  due  to  the  velocity  through 
X 

the  water  of  Ap,  as  shown  in  Figure  Dl. 
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The  Incremental  force  on  an  element  of  length  Ap  Is: 

AF.  »  K(u.  +  0).  x  o)3Ao 
A  A  A 

The  Integral  of  the  incremental  forces  is  as  follows: 

V2 

<F«cc>X  ‘  K  I  R  +  2“X  '  +<"x  *  F)3]  dp 

V2 

The  above  integral  may  be  simplified  if  we  observe  that: 

2ux  *  ^“x  x  p)  *  2\  '  x  s^p)  m  f(ux»u,x)p 

(<«x  x  p)J  -  u^p3sin(p,u>x)  =  ui^p3, 

since  by  hypothesis  sin  (p,u>  )  *«  ±  i.  Here  u;  is  taken  to  be 

K  K 

normal  to  the  vector  p  . 


vector 
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Also  noting  that 

V2  +V2 

j  pdp  *  0  ,  and  J  p2dp  ■  d^/12  » 


'V2 


"V2 


then  if  (Fftp  be  the  approximate  force  obtained  by  neglecting 
the  angular  velocity  ou  ,  we  have 

A 

+Y2 

'Vr1  I  i  dp"Kdx 

-dx/2 

The  relative  error  is 


(f.pp'x 


a/12)  c^y^r1 


so  that  if  the  velocity  Oue  to  rotation  or  one  end  of  tne  element  with 
respect  to  the  other  were  equal  to  the  velocity  of  the  midpoint  through 
the  water,  the  emission  of  rotation  causes  an  error  of  8.25  per  cent. 

If  the  relative  velocity  due  to  rotation  were  one-half  of  the  translation 
velocity,  the  error  would  be  about  2  per  cent.  This  can  always  be 
achieved  by  using  short  cable  elements,  but  apparently  this  is  not  a 
stringent  limitation.  Therefore  when  calculating  the  velocity  forces,  the 
rotation  of  the  element  will  be  ignored. 


By  a  direct  three  dimension  generalization  of  the  normal  and  tan¬ 
gential  forces  discussed  in  Reference  2,  one  can  obtain  a  representation 
of  the  tangential  force  with  a  magnitude  of  d^  and  direction's^; 

and  a  representation  of  the  normal  force  with  a  magnitude  of 

d*<K  *ln\  *  (-“x  x  VM>3 "  VS\ x  V 

with  a  direction  of  q^  . 

The  moment  above  the  c.m.  due  to  lift  is  now  computed.  Let  v  be 

the  velocity  of  Ap  through  the  water,  then  v  =  u  +  uj  x  p's  .  The 

A.  A  A.  , 

incremented  force  experienced  by  the  element  becomes: 


AF  *  QsS?  sin2  (v,-s  )Ap 
A, 

where  v  x  (-s  )  *  v  sin  (v,-s  )  .  Then  AF  takes  the  form 

K  \ 

AF  =  Cs  | s  x  v|2  Ap  .  The  element  moment  about  the  c.m.  of  the  cable 
element  takes  the  form: 

j  AM  -  ps^  x  (AF)qx  -  p  Qa  |  s^  x  v(2  x  Ap 

I  since  x  ?  «  s2v2  -  (s^  .  v)2  ,  then  an  integration  from  -d^/2 

to  d  /2  yields: 

* 

|  V2 

\“4(VV  I  [p^  -  («x  .  v)3]dp 
-<i/2 

Noting  that  (s^  .  v)2  -  (s^  .  u^)2  ,  and  also  that  integration  of  terms 
j  with  odd  powers  of  p  will  give  no  net  contributions,  the  above  expression 

becomes: 

i  _  _  V* 

V  ■  x  \  I 

■\n 

+V2 

■  ^’x  x  ’x(“x  ■  \  *  *x>  J  2p3<ip 

•V2 

Finally  the  moment  of  the  lift  about  the  c.m.  reduces  to  the 
form 

V  ■  «K/6>  <“x  •  \  *  *x>“x  *  \ 
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APPENDIX  E 


ACCELERATION  FORCES 


In  this  appendix  there  are  two  mutually  exclusive  but  alternative 
physical  assumptions  to  be  made  regarding  the  composition  of  these 
forces.  The  first,  which  we  label  Approach  A,  is  a  three-dimensional 
formulation  of  the  physical  assumption  embodied  in  Reference  1.  The 
second,  which  we  label  Approach  B,  is  a  Kirchhoff's  general  formulation 
of  the  added  masses  of  a  submerged  body. 

The  mathematical  expression  of  Approach  A  is  simpler  than  that  of 
Approach  B.  Approach  A  is  always  used  for  the  cable.  Approach  B  may  be 
more  desirable  for  the  vehicle.  Equations  are  derived  for  each  approach. 

These  forces  and  moments  contain  the  accelerations.  The  accelera¬ 
tions  are  the  unknowns  for  which  the  final  linear  equations  are  to  be 
solved.  Approach  A  and  Approach  B  assign  different  values  to  the  co¬ 
efficients  of  these  equations,  and  therefore  differ  basically. 


APPROACH  A 

In  Approach  A  we  are  concerned  only  with  the  force  due  to  the  in¬ 
duced  acceleration  of  the  fluid.  As  in  Reference  1  we  assume  that  this 
force  is  proportional  to  the  component  of  the  acceleration  normal  to  the 
cable  element,  and  that  there  is  no  moment  due  to  the  acceleration. 


FIGURE  El.  UNIT  VECTORS  s  ,  q  ,  AND  ACCELERATION  a 

XX  X 


Let  n^  be  the  unit  vector  normal  to  and  in  the  plane  of  a^ 

and  s  ,  as  shown  in  Figure  El.  The  acceleration  a  is 
X  X 

•x  ■  C\  ■  \\  *  C\  ■ 
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Then  the  component  of  acceleration  normal  to  the  cable  element  is: 

c\  ■  \>\ m  \  '  <\  :  V'x  ■ 

Hence  the  force  due  to  acceleration  of  the  element  is 

\  •  "x  ■  •e’m<‘x  •  *x  •  Vx> 

‘  -ex<*x  -  *x  •  *x  V 

where  e'm^  -  e^  is  a  parameter  which  depends  upon  the  added  mass  of 
fluid  of  cable  element. 

APPROACH  B 

Consider  a  body  moving  in  a  perfect  fluid  which  is  otherwise  at 
rest.  The  motion  of  the  fluid  is  assumed  to  be  irrotational.  The  hydro- 
dynamic  forces  and  moments  experienced  by  this  body  due  to  its  linear  and 
angular  accelerations  are  given  by  the  following  (Reference  3): 


-"[  l 

j-i 


l 

3*1 


]• 


(1-1,2, 3) 


Mi- 


j-1 


i+3,J 


3 


(t)n 


i+3 ,  j+3 


(i  -  1,2,3) 


j-1 


u”(t)  are  the  components  of  the  translational  velocity  of  the  origin  of 

a  system  of  rectangular  axes  fixed  to  the  body,  and  u)j  are  the  components 

of  angular  velocity  of  the  body  about  this  origin  as  shown  in  Figure  E2. 
The  density  of  the  fluid  is  given  by  p,  and  p  are  the  "added  masses". 

The  conctants  p  are  determined  exclusively  by  the  shape  and  position 

of  the  body  surface  relative  to  the  moving  coordinate  axes  fixed  to  the 
body.  It  is  also  known  that  *  p^.  consequently  the  above  equations 

for  forces  and  moments  involve  21  possible  independent  constants 
(Reference  3).  If  the  surface  of  the  body  has  one  or  more  planes  of 
symmetry,  the  number  of  independent  constants  a*e  reduced.  For  example, 
if  the  -  plane  is  a  plane  of  symmetry,  it  can  be  proven  that 

21  independent  constants  reduce  to  12.  If  the  53?3  -  plane  is  also  a 
plane  of  symmetry,  then  the  12  independent  constants  reduce  to  8.  If 
there  is  also  symmetry  of  body  surface  about  the  £1^3  -  plane,  the 


-56- 


8  independent  constants  reduce  to  6,  so  that  only  the  constants 
U^^(i  ■  1,  . ..,6)  remain.  However,  In  this  analysis  ve  chall  retain 

all  the  constants  or  "added  masses." 


*3 

FIGURE  E2.  FIXED  FRAME  ,  AND  BODY  FIXED  FRAME  5^  5, 


Let  5i »  5b  >  5s  denote  the  components  of  linsar  acceleration 
of  the  moving  origin  in  the  direction  of  the  axes  fixed  to  the  body. 

Also  let  F-  and  M_  denote  the  components  of  hydrodynamic  forces  and 

*i  h 

momsnts  respectively  in  the  direction  of  the  axes  fixed  to  the  body. 

The  above  equations  for  the  forces  and  moments  now  take  the  following 
form: 

\  ’  -  •>  Lclij«j  -  C3iAjl 

\  ‘  *  11  K‘j  -  C2u“!)] 

where  Cl ^  ^  for  i  s  3,  J  i  3;  02^  -  H1+3j+3»  for  i  +  3  >  3, 

j  +  3  >  3;  03^  ■  j+j,  for  i  S  3,  J  +  3  >  3,  in  the  equation 
for  the  force;  and  C3^  »  i  +  3>3,  J  S3,  in  the  equation 

for  the  mossint. 
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To  refer  the  forces  and  moments  relative  to  the  directions  of  the 
inertial  axes  we  note  that  for  any  vector  we  can  write  ■  c^  Vg  , 

i.e. t  Vt  «  ctj  V’ ,  and  -  c^  Vk<  Here  the  prime  Indicates  evaluation 
with  respect  to  the  moving  frame,  as  shown  in  Figure  E3. 


FIGURE  E3.  UNIT  VECTORS  ON  FIXED  AND  MOVING  FRAME 


Referring  V  to  the  basis  i,  J,  k  and  V'to  the  basis  l,  m,  then 
dWdt  -  ^V'/dt  +  x  V.  In  particular  du>/dt  »  diu'/dt  +  (5  x  u>  -  is'/dt, 
so  that  a)  -  c.  .u),  . 

*i  J  Sj 

Hence 

\  ‘  CIJ  %  '  *  >  ciJ  [Cl«5k  -  c3ik\J 

"  "  p  clj  [C1jk5k  "  C3jkciklLxi] 

•  • 

where  u>g  ■  Cjk»»x  >  *nd 

\  ■  ’  cii  [cVk  -  “jk'iA,] 

The  next  step  in  the  analysis  is  concerned  with  the  expression  for 
§k  in  terms  of  the  inertial  coordinates  x1  ,  % ,  Xa .  From  Figure  E3,  and 
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also  from  a  coordinate  transformation  we  write  =  c^j5j  +  where 

d:  ,  dg ,  dg  are  the  ^ ,  Xg  ,  Xg  coordinates  of  the  origin  of  the  moving 
frame.  The  transformation  is  also  orthogonal  so  that  c  c^  *  6^* 

c..c.  .  =  6,,  ,  and  5.  *  c. ,(x  -  d,).  We  also  note  that 
ij  kj  ik’  j  ij  1  i 

r  =  xl  +  yj  +  zlc  =  +  T]m  +  Jn  , 

and  then  by  the  usual  rules  of  kinematics  we  obtain: 


and 


r  *  xi  +  yj  +  zk  *  +  T]m  +  £n  +  to  x  r, 

1 1  • »  • 
r  «  xi  +  yj  +  zk 


■  §i  +  Tjm  +  £n  +  2u>  x  (§4  +  Tim  +  £n)  +  coxrxcox(<Dxr) 
This  equation  can  be  rearranged  in  the  following  form: 

§  -  ?4  +  T{tn  +  Cl)  *  r  -  2u)  x  r  -  (u)  x  r)  +  uu  x  (u>  x  r) 

-  AA  +  AB 

. -  JU  JL  —  ••  • 

Where  AA  »  r  •  ui  x  r.  so  that  AA,  *  x  -  e.  u)  x  ; 

J  j  jpq  p  q 

and  15  -  -2u)  x  r  +  <0  x  (co  x  r) ,  so  that  AB^  *•  -  ejpq2oUpXq 

U 

+  e,  uu  e  .  U)  x„.  Thus  the  k  th  component  of  §  is: 
jpq  p  qst  s  t 

\ m  V"j +  *V 

•  •  •  • 

*  c..  (x  -  e.  0)  x  )  +  c,,  (e,  e„.u> ji)  x.  -  2c,  cu  x  ) 

Ik'  1  ipq  p  a'  Ik'  jpq  qst  p  8  t  Jpq  P  q7 


where  the  first  term  on  the  right  hand  fide  contains  accelerations,  the 
second  term  does  not  contain  accelerations,  and  x^  is  the  q  th  component 

of  r  where  r  is  the  radius  vector  from  the  inertial  origin  to  the  origin 

of  the  moving  frame.  Finally 

•  •  ♦  • 

F  =  -p  c,,  Cl.,  c  .  (x  v-e — co  )  -  p  c,,C3,.c..u) 

xt  H  ij  jk  rk'  r  rpq  p  q'  ij  jk  4k  x£ 

-p  c.  ,  Cl,,  c  ,  (e,  e  co  co  x„  -  2e  co  x  ) 
y  ij  jk  rk'  jpq  qst  p  s  t  rpq  p  q' 

•  *  «  • 

H  =  -  p  c,,C3,,  c  ,(x-e  u>x)-p  c  ,C2  ,  c  ,.co  ) 
xt  y  ij  jk  rk'  r  rpq  p  q'  ij  jk  ik  x^ 

-pc,  ,C3,,  c  ,  (e,  e  cu  uj  xt  -  2e  u»  x  T 

p  ij  jk  rk'  jpq  qst  p  s  t  rpq  p  q 

(Reverse  Page  60  Blank) 
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APPENDIX  F 


FIN  FORCES  AND  MOMENTS 


The  vehicle  considered  here  has  a  pair  of  fins  in  a  vertical  plane 
(one  fin  above  the  vehicle  and  one  below),  and  a  pair  of  fins  in  the 
horizontal  plane  (one  to  port  and  one  to  starboard).  These  fins  are 
subsequently  referred  to  as  the  vertical  fins  and  the  horizontal  fins, 
respectively. 


Consider  first  the  vertical  fins.  Figure  FI  illustrates  the  forward 
half  of  one  of  the  vertical  fins.  The  velocity  through  the  fluid  of  the 
vertical  fins  (see  Appendix  A)  is: 


vf 


xv+l  +  (dv+Vdv+P 


v+2 


•  XH-1) 


3 

FIGURE  FI.  VERTICAL  FIN  AND  ITS  VELOCITY 


The  lift  (or  sidewise  force  in  the_case_o£  vertical  fins)  will  be 
*  determined  by  the  general  .relation  I  »  T  x  pU,  where  L  is  the  lift, 

F  is  the  circulation,  and  U  is  the  velocity  of  the  fins.  The  vector 
F  will  be  assumed  to  be  equal  to: 
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(constant) (velocity) (sine  of  the  angle  of  attack) (unit  vector) .  The 
effect  of  spanwise  flow  of  the  fluid  will  be  ignored  in  this  formu¬ 
lation.  This  assumption  implies:  (1)  that  the  circulation  vector  is 
parallel  to  the  axis  of  the  fin,  i.e.  ±  n;  (2)  that  the  angle  of 
attack  is  measured  in  the  §i§2  plane,  as  shown  in  Figure  F2;  and 
(3)  that  the  velocity  with  which  we  are  concerned  is  v^  -  (v^  .  n)n  , 

i.e.  the  projection  of  v  ^  on  the  §2  plane. 


vvf  ’  (vvf  *  n)n 


FIGURE  F2.  FIN  AND  VELOCITY  IN  Si§a  PLANE 


The  angle  of  attack  is  a  +  (-c^),  where  a  is  the  kflown _de f l£c t ign 

of  the  fin  relative  to  the  vehicle,  where  or,  =  arctan  (m  .  v  ,  /l  .  v  , 
_  _  _  _  1  vfp  _  vfp 

arctan  (m  .  v^/X  .  v  ^).  Hence  the  circulation  vector  becomes  T  * 
(constant)  v^n  3in  (a  -  c*i )  and  the  aJ,  dewise  force  takes  the  form: 


(constant)vyfp  sin(of  -  Qi)n  x  p(vyf  -  (vyf  .  n)n) 

where  v  is  the  magnitude  of  v  ,  .  Let  q  , „  uc  a  unit  vector  in  the 

vfp  vfp  ^v+2 

direction  vyj,  and  q,^  another  unit  vector  defined  by  n  x  . 

Since  n  x  n  =  0,  the  vector  factor  in  the  above  equation  n  x  p  v  ^  * 

pvvfqv+38irKn,Vvf)  *  p  ^<1^3  •  Then  the  lift  *  (constant) 

pv2^  q^^  sin(a  *  Oi ) .  Since  the  fin  is  of  low  aspect  ratio  and  is 

in  the  neighborhood  of  a  large  body,  the  lift  may  be  written  in  the 
following  form: 
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FLVF(o  "  04  ^  'vfp  **v+3 

where  Cj  is  an  empirical  constant  greater  than  zero,  and  ^VF ~ 

is  an  empirical  odd  function  of  its  argument,  the  angle  of  attack. 

Both  and  may  be  obtained  experimentally. 

If  a  is  positive  when  the  forward  edge  of  the  fin  is  displaced 
to  starboard^  and  F^  is  positive  when  its  argument  is  positive,  then 

positive  lift  (sidewise  force)  tends  to  force  the  vehicle  to  starboard 
and  to  increase  the  ^  coordinate.  The  drag  on  the  vehicle  13 

similarly 

-C*  Fcvf  (or  -  Qi>  <fpV2 

where,  however,  is  an  even  function  of  its  argument.  The  moment 

about  the  center  of  gravity  of  the  vehicle  due  to  the  vertical  fins  is: 

*Vvfp^  FLVF^®  "  Q^C»(dv+4  "  dv+2^*vfl  X  qv+2 

^vvfp^*  C*  *v+l  x  qv+2  FDVF  ^  dv+2^ 

which  reduces  to  the  following  form: 

<vvfp)3  (dv44  “  dv+2^*v+l  x  ^  V3  FLVF(a  '  01 ) 

**  qv+2  FDVF^®  "  ^ 

By  proceeding  in  an  analogous  manner,  similar  results  are  obtained 
for  the  horizontal  fins.  In  this  cas»  the  velocity  of  the  horizontal 
fins  is: 

vhf  "  xv+l  +  (dv+3/dv+l'(xv+2  ‘  Xv+P 

Let  q^^  be  a  unit  vector  in  the  direction  of  vvf;  *  -  m  x  q^; 

vhfp  t*ie  ProJect^on  °*  on  c^e  5i5s -Plane;  Pj  ■ 
arctan  (n  .  v^^/i  .  v^) ;  and  let  P  be  positive  when  the  leading  edge 
of  the  fin  is  deflected  down. 

The  lift  due  to  the  horizontal  fins  is  then 

**  *vhfp^  Vf5  FLHF^  * 
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where  Cb  and  Fjjjp(P  "  Pi )  are  analogous  to  Cj  and  F^,  respectively. 
Therefore  P  >  0  tends  to  increase  3.  The  drag  on  the  horizontal 

fins  is  _ce  (vhfp^3qv+4  FDHF^  "  ) »  and  the  moment  about  the  center 

of  gravity  of  these  forces  is: 

(vhfp^dv+3  '  dvfPSv+l  X  ^  qv+5  FLHF(P  &■) 

^  qv+4  FDHF^ 

The  total  force  FIN,  and  moment  FINM  acting  on  the  vehicle  are  the 
sums  of  the  above  forces  and  the  sums  of  the  above  moments.  They  are: 

FIN  -  (vvfp)8[C9  V3  Flvf  (a  -  Oi)  -  CA  FDVp(of  -  as,)] 

a.  'v2  r n  v  to.  -  o  %  -  r  v  (a  -  a  \1 

'  ’hfp'  ‘‘v+5  "LHF  'r  "A/  ■‘v+4  '  DHF'r  rl  7  J  ’ 

and 

™  *  (\fp)2(dv44  ’  <W%+1  X  CQj  qv+3  FLVF(a  ’ 

C4  qv+2  FDVF^a  ”  01  ^ 

+  (vhfp)2 (dv4-3  "  dv+2)sv+l  X  ^  qv+5  FLHF  '  Pl ) 

^v+4  FDHF  ^  ^  ^ 
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APPENDIX  G 


COORDINATE  TRANSFORMATION  BETWEEN  INERTIAL 
COORDINATES  AND  COORDINATES  FIXED 
IN  THE  VEHICLE  WITH  ORIGIN  AT 
THE  CENTER  OF  GRAVITY 


The  coordinates  fixed  in  the  vehicle  with  their  origin  at  c.g.  of 
vehicle  are  denoted  by  §,  1),  £.  The  coordinates  x,y,z  denoted  the 
coordinates  in  a  frame  which  is  parallel  to  the  inertial  frame.  The 
origin  of  both  frames  is  at  0,  and  the  £,  T]  plane  contains  the  hori¬ 
zontal  fins,  as  shown  in  Figure  Gl.  This  appendix  presides  the  co¬ 
efficients  c^j  to  transform  vector  components  between  the  two  frames. 


FIGURE  Gl.  INERTIAL  AND  BODY-FIXED  COORDINATES 


For  a  linear  orthogonal  transformation  x^  1 
the  coefficients  c^  must  satisfy  the  following 


'  cijSj*  (®ee  Appendix  B) 
orthonormal  relations: 


CijCik  *  6 jk  »  cjicki  *  6kj 
The  transformation  in  terms  of  matrices  is: 


-65- 


If  §  =  1,  1)  =  0,  £  =  0,  are  the  components  of  -s,  then  the  substitution 
of  these  coordinates  in  the  preceding  matrix  transformation  gives 


X  = 

cn 

.  (1.0)  =  -s 

X 

y  = 

C21 

00 

II 

o 

H 

N— f 

z  = 

C31 

.  (1.0)  =  -sz 

If  next  we  let  ?  =0,  T|  =  1, 

C  =  0, 

a  substitution  in  the  transformation 

with  condition  of  zero  roll  gives  z  =  c^.Cl.O)  *  0. 

With  the  above  mentioned  quantities  for  x,y,  and  z,  the  coefficient 
transformation  takes  the  following  form: 


the 

c33 


The  crthonormal  relations  are  applied  to 
c|jS  and  the  s^s.  For  example,  * 


find  the  relation  between 
1  ■  s3  +  Cj3  .  Solving  tor 


we  get  c233 


s2  +  s2  ,  so  that  c_ 
x  y  33 


when 


s^  =  0,  z  =  c33C>  80  that  we  take  the  positive  square  root. 


'33 


Js2"+ 


1/X 


where  >  is  defined  for  convenience  to  be  1/ 


•Js2  +  s3  . 


hence 


For  the  second  and  third  rows  the  orthonormal  relations  give 

s  s  +  c«,  c0,  =  0,  thus  c„.  =  -s  s  /c-~  =  -s  s  X  . 
y  z  23  33  *  23  y  z  33  y  z 

For  the  first  and  third  row  the  orthonormal  relations  give 

Vz  +  C13C33  =  °’  hence  c13  =  -Sx8z/c33  =  “VzX  ‘ 

Also  c22  +  C22  =  1,  and  -s^^  "  syc22  =  0»  are  consequences  of 
the  orthonormal  relations.  The  last  condition  yields  s2  c2  =  s3  c3 


Using  this  relation  and  the  other  condition  gives: 


o2  -  o2^2  jl  r>2  _  _2  2  ,  2  2 

s  -  s  c10  +  s  c00  =  s  co0  +  s  c00 
x  x  12  x  22  y  22  x  22 

from  which  ct0  =  s2/(s2  +  s2),  hence  c00  =  u\s  ,  where  u  denotes  +  1  or 
22  x  x  y  22  x 

-  1.  In  addition 

C12  =  (-Sy/Sx)C22  *  -UX8y 
The  determinant  of  the  transformation  is 


-s 

X 

-^X8y 

-s  s  ) 
X  z 

-s 

y 

W8x 

-s  s  X 
y  z 

-s 

7 

0 

(1/X) 

which  equals  -ias2  -^(s2  +  8y' •  Since  the  determinant  must  be  equal  to 


s2  -pi(sa  -  -2' 

Z  X 

+1,  then  it  follows  that  ^  ■  -1. 

In  summary,  Xj^  -  cij^j’  *  cijxi  »  where 
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